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ABSTRACT 


In  an  acoustic-oceanographic  experiment  (AFAR)  performed  near  the  Azores  in 
1975,  frequencies  from  400  to  4670  Hz  were  transmitted  over  a  35-km  wholly  refracted 
path.  In  addition,  a  separate  data  set  was  gathered  over  a  3-km  path  in  1973.  Meas¬ 
urements  of  the  intensity  coherence  function  of  time  and  frequency,  intensity  moments, 
and  the  probability  distribution  of  intensity  are  presented  and  compared  with  predic¬ 
tions  of  fluctuations  due  to  internal  waves. 
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INTRODUCTION 

In  this  paper  we  compare  with  theory  some  observations  made  by  Ellinthorpe1-2  of 
the  statistics  of  acoustic  intensity  observed  at  the  Azores  Fixed  Acoustic  Range 
(AFAR).  We  treat  two  separate  experiments,  one  in  1973  and  a  second  in  1975.  Our 
theoretical  predictions  for  the  intensity  statistics  are  obtained  from  the  path-integral 
formulation  under  the  assumption  that  the  effects  of  internal  waves  dominate  the 
acoustic  fluctuations.  This  approach  is  presented  in  Dashen3  and  Flatte'  et  al,4  and 
reviewed  in  Flatte'.5  A  companion  paper6  presents  the  pertinent  theoretical  considera¬ 
tions  relevant  to  the  comparison  between  the  predictions  and  observations  of  intensity 
at  AFAR.  In  an  earlier  pair  of  papers  7  8  a  similar  comparison  for  the  second  moment 
(mutual  coherence  function)  and  a  detailed  description  of  environmental  observations 
were  given.  Here,  we  present  the  fourth  moments  (more  specifically  <I(0)I(At)>  and 
<I(0)I(A<r)>;  where  A t  and  A a  are  separations  in  time  and  frequency,  and  the  scintilla¬ 
tion  indices);  higher  moments;  and  probability  distributions. 

It  is  now  accepted  that  ocean  acoustic  fluctuations  having  time  scales  from 
minutes  to  days  are  caused  predominantly  by  internal-wave-induced  oceanic  variability 
in  the  sound-speed  field.5  Ocean  experiments  have  sampled  the  fluctuating  acoustic 
intensity  field  at  a  variety  of  locations  and  ranges.  Adopting  the  A—  $  scattering 
parametrization,  we  note  three  regions  that  catagorize  acoustic  fluctuations  (Figure  1). 
In  the  unsaturated  regime  (A^2  <1  or  $<1),  where  the  intensity  field  is  expected  to 
obey  log-normal  statistics,  there  are  no  experiments.  Both  the  San  Clemente  experi¬ 
ment,9  and  the  experiments  at  Cobb  Seamount  10  11  are  near  the  A  Or  =1  boundary.  All 
these  experiments  found  stronger  intensity  fluctuations  than  one  would  expect  over  a 
weakly  scattered,  unsaturated  path.  In  the  saturated  regime  (A<&>1)  where  the  com¬ 
plex  field  components  are  expected  to  be  Gaussian,  there  are  the  MIMI  experi¬ 
ments,41213  and  an  experiment  southwest  of  Bermuda  14  The  ENACEX  observations15 
sampled  three  ranges  with  a  frequency  of  460  Hz  and  several  paths  and  thus  sampled  a 
broad  range  in  A$  -space.  ENACEX  was  used  to  estimate  several  fluctuation  parame¬ 
ters  directly  from  the  measurements  The  paths  predicted  to  be  nearly  saturated  exhi¬ 
bited  partially  saturated  behavior  Since  the  time  of  these  experiments,  the  path- 
integral  theory  has  developed  sufficiently  so  that  intensity  observations  where  the 
environment  is  adequately  known  can  be  used  to  test  theoretical  predictions  with 
greater  precision 

The  AFAR  propagation  occurred  over  two  paths,  one  of  nearly  3  km  in  1973  and  a 
fixed  path  of  35  km  during  1975  These  observations,  for  the  several  center  frequencies 
transmitted  at  each  range,  place  the  short  range  experiment  in  the  unsaturated  regime 
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and  the  35  km  path  in  partial  saturation  (Figure  1).  Simultaneously  with  the  1975 
experiment,  an  extensive  set  of  environmental  measurements  were  made. 

The  paper  is  divided  into  three  parts.  Part  I  reviews,  from  our  earlier  paper,  the 
acoustic  experiments  and  the  oceanographic  observations  determining  the  parameters  of 
the  internal-wave  model  (we  use  the  GM-formalism).11  We  also  introduce  the  theoreti¬ 
cal  framework  discussed  in  our  companion  paper.6  The  intensity  measurements  and  the 
fourth  moments  are  presented  in  Part  II  and  the  probability  densities  and  higher 
moments  in  Part  III.  The  paper  concludes  with  a  summary:  Part  IV. 
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I.  EXPERIMENT  REVIEW  AND  OVERVIEW  OF  THEORY 

i 

A  brief  overview  of  the  experiment  is  presented  here.  The  interested  reader  may 
find  more  experimental  detail  and  an  associated  list  of  references  in  our  earlier  paper.8 

|  Several  different  frequencies  were  transmitted  over  each  of  two  refracted  paths  from  a 

fixed  transmitter  at  540-m  depth  located  west  of  Santa  Maria  Island  in  the  Azores.  A 
2798-m  path  (inclined  upward  at  5  degrees)  was  monitored  for  seven  hours  in  November 
1973;  the  receiver  was  a  hydrophone  suspended  at  325-m  depth  beneath  a  trimooring. 

i  Three  pulses  having  center  frequencies  749.8,  1499.8  and  2800.0  Hz  were  transmitted 

t 

j  every  1.048576  s.  A  longer  35026-m  path  was  monitored  with  a  fixed  receiver  located  at 

736-m  depth  for  over  186  h  beginning  in  March  1975,  simultaneously  with  an  extensive 
oceanographic  experiment.  Three  130  ms-pulses  of  412.46,  1012.33  and  4673.01  Hz,  and 
a  fourth  1.25-ms  3199. 5-Hz  pulse,  all  having  Gaussian  envelopes,  were  transmitted  at  a 
10.15808-s  repetition  interval.  In  the  language  of  communication  theory,  the  short 
pulse  is  used  to  measure  the  channel  frequency  response  and  the  long  pulse  measures 

the  time  response. 

i 

|  Environmental  measurements  obtained  during  the  1975  experiment  determine  the 

background  sound-speed  field  and  the  parameters  of  the  Garrett-Munk  (GM)  internal- 
wave  model.16  The  sound-speed  profile  has  a  minimum  near  714  m  in  addition  to  the 
deeper  SOFAR  axis.  The  intervening  maximum  near  850  m,  is  caused  by  the  intrusion 

1 

j  of  Mediterranean  water.  The  714-m  minimum  creates  a  duct  for  the  long-range  ray, 

[  which  has  two  turning  points  and  two  caustics.8  Specification  of  the  internal-wave  field 

j  involves  parametrization  of  the  rms  vertical  displacement  bandwidth  j,,  and  the 

scaling  parameter  n0B  of  the  dispersion  relation.  Using  a  combination  of  horizontal, 
vertical,  and  temporal  environmental  data,  values  for  the  parameters  were  obtained  8 
of  n0B  =3  0  cph-km,  ^  =13.1  m  and  j *  =3.  These  are  the  same  as  the  canonical  GM 
model  with  the  important  exception  that  ^j,  the  parametrization  of  the  internal-wave 
energy  level,  is  much  larger  than  the  canonical  value  of  nearly  7  m  and  is  relatively 
high  compared  with  other  oceanic  internal-wave  observations.17  For  our  purposes  this 
means  AFAR  is  a  good  place  to  study  internal-wave-induced  acoustic  fluctuations. 

With  the  model  sound-speed  statistics  specified,  then  the  theoretical  predictions  for 
the  intensity  statistics  can  be  obtained  for  the  AFAR  site.  The  companion  paper 
addresses  the  derivation  of  the  expressions  using  the  path-integral  formalism  6  The 
theoretical  quantities  of  interest  include  the  intensity  correlation  function  in  time 
K(At)  and  the  micropath  bandwidth  function  from  which  the  intensity  correla¬ 

tion  of  frequency  is  obtained.  For  these  quantities  the  comparison  between  observations 
and  theory  is  contained  in  the  following  intensity  correlations. 
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where  the  brackets  <>  indicate  an  average  over  the  duration  of  the  experiment,  which 
is  taken  to  represent  an  ensemble  average  over  the  statistical  internal-wave  field.  (  The 
mean  intensity  is  by  definition  equal  to  unity.)  The  solution  of  K(Af)  involves  numerical 
integration  of  a  second-order  differential  equation  along  the  equilibrium  ray,  as 
described  in  the  companion  paper.6  The  theory  for  Q  was  described  earlier.7 

In  full  saturation,  K  (A<)  and  for  small  time  lags,  D(At)  «  (i/At)2.4  In 

an  analogous  manner,  in  the  partially  saturated  regime,  an  expression  for  K  can  be 
defined  for  small  time  lags, 

K(A0  (3) 

This  defines  an  intensity  decorrelation  time  scale  l/v',  which  can  be  calculated  by 
means  of  ordinary  integrals  along  the  ray  from  source  to  receiver,  for  given  internal- 
wave  parameters.8 

The  behavior  of  |Q(A  cr)|2  can  be  approximated  for  some  purposes  by 

|Q(Aa)|2«  e  ~(r<^’2  (4) 


where  r0  can  be  calculated  by  means  of  ordinary  integrals  along  the  ray  as  well,  again  if 
the  internal-wave  field  is  quantitatively  characterized.8 

For  the  moments  of  intensity,  as  shown  by  Dashen18  the  general  asymptotic  form 
is: 


<1*  >  =  /VJexp 


(5) 


The  form  of  (5)  for  small  q  with  /2(q)  =  1  is  used  here: 


<lv  >  =  .V!  exp 


MAM) 

i  1 


(6) 


where  the  value  for  q  is  obtained  from  the  effects  of  internal  waves  by  the  path-integral 
technique.19  For  small  q,  (6)  takes  on  the  more  familiar  form 


<IV  >  =  .V 


1  + 


,V(,V-t) 


(-) 
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The  moments  obtained  from  the  K-distribution20  are  also  considered,  where  we  use  the 
form 

<I^>  =  ^!lLX.+  l/7hN,  (8) 

Not  only  the  solutions  K(Af)  and  Q(A t)  but  also  all  other  quantities  including 
v,  v tq  and  7  are  obtained  by  numerical  integrations  along  the  deterministic,  back¬ 
ground  ray  from  source  to  receiver.  The  AFAR  background  ray  path  is  calculated  using 
a  sound-speed  profile  obtained  from  an  average  of  65  STDSV  casts  taken  during  the 
1975  field  experiment.  There  is  an  issue  of  determining  the  mean  background  profile 
and  distinguishing  it  from  the  fluctuating  profile.  Even  with  the  formation  of  an  average 
profile  from  many  casts,  the  intensity  related  quantities  including  K,  A,  and  Q  are  func¬ 
tions  of  how  the  background  profile  is  smoothed.  The  smoothing  (and  interpolation) 
algorithm  used  here  is  a  cubic  spline  technique,21  involving  a  parameter  which  is  the 
rms  deviation  of  the  smooth  profile  from  the  measured  points.  The  predictions  are  sensi¬ 
tive  to  this  parameter,  as  was  described  in  our  earlier  paper8  involving  Q.  For  the  com¬ 
parison  with  AFAR  the  smoothing  parameter  is,  unless  specified,  set  to  1  cm/s.  This 
level  of  smoothing  effectively  filters  out  much  of  the  high  vertical  wavenumber  variabil¬ 
ity  (smaller  than  20  m)  from  the  background  sound-speed  profile. 


V 


n.  FOURTH  MOMENTS 

In  treating  the  statistics  of  any  time  series  it  is  useful  to  initially  examine  the 
character  of  the  series.  This  is  helpful  at  AFAR  where  there  are  intensity  measurements 
from  experiments  in  two  quite  different  scattering  regimes.  Figure  2  shows  time  series 
from  three  frequencies  at  each  of  the  two  ranges.  The  entire  available  data  set  is  shown 
for  the  2.8-km  path;  20  h  of  data  (an  inertial  cycle  at  the  AFAR  latitude)  is  shown  from 
the  35-km  path.  Note  that  the  short-path  intensity  ordinate  is  a  factor  of  2  smaller 
than  that  for  the  longer  path. 

Over  the  shorter  path,  the  data  series  from  the  different  frequencies  display  the 
same  general  oscillation  sequence  in  time.  .As  the  acoustic  frequency  increases,  the  oscil¬ 
lations  increase  in  amplitude.  This  is  the  behavior  expected  in  a  weak,  geometric 
regime,  where  theoretically  the  scintillation  in  log-intensity  scales  with  frequency.  The 
observed  scintillation  indices  for  log-intensity  (and  intensity)  are  0.10  (0.10),  0.21  (0.15), 
and  0  59  (0.34)  for  750,  1500,  and  2800  Hz  respectively. 


'i 


£ 


i 

* 
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Over  the  35  km  path,  at  longer  time  scales,  the  data  behave  in  a  similar  manner. 
For  example  there  is  a  broad  increase  in  intensity  centered  about  10  hours  and  again 
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after  16  hours.  However,  comparing  the  shorter-time-scale  events  oscillating  inside  the 
long  period  envelopes,  the  three  series  are  quite  different  and  the  high  intensity  spikes 
between  frequencies  are  often  not  in  alignment.  Over  very  short  time  scales  the  series 
for  the  different  frequencies  are  quite  different.  The  observed  scintillation  indices  of 
intensity  for  the  35  km  path  are  1.24,  1.21,  and  1.47,  for  410,  1010,  and  4670  Hz  respec¬ 
tively.  The  corresponding  log-intensity  variances  are  1.8,  1.7,  and  1.9.  These  are  above 
the  saturation  limit  of  1  for  intensity  and  1.66  (5.6  dB)  for  log-intensity.  This  is  the 
behavior  characteristic  of  partial  saturation,  in  which  the  deterministic  ray  has  been 
broken  apart  into  many  microrays  by  the  internal-wave  fluctuations.  The  level  of  scin¬ 
tillation  at  4670  Hz  is  higher  than  expected  on  a  curve  asymptotically  approaching  1 
from  above.  (The  approach  to  1  is  very  gradual  in  partial  saturation.) 

We  turn  now  to  the  comparison  of  the  observed  fourth  moments  (intensity  correla¬ 
tions)  and  theory.  A  quantitative  examination  is  made  of  the  observations  over  the  35- 
km  path.  Because  of  mooring- motion  problems,  the  2.8  km  path  cannot  be  used  with 
the  intensity  correlation  studies.  The  temporal  correlation  function  r(Al)  for  1010  Hz 
is  plotted  in  Figure  3.  Standard  FFT  techniques22  for  the  correlation  functions  are  used 
and  the  confidence  limits  are  obtained  from  the  equations  in  Appendix  B  of  Spiesberger 
and  Worcester.14 

The  prediction  K(Al)  is  a  function  of  smoothing  of  the  background  sound-speed 
profile.  K(Af)  curves  for  two  smoothings,  1  and  1.5  cm/s,  are  shown  in  Figure  3;  90% 
confidence  intervals  are  shown  for  a  few  lags  on  the  dotted  data  curve.  As  one  might 
expect,  the  greater  the  smoothing,  the  more  slowly  the  predicted  correlation  function 
falls  with  lag.  There  is  a  smooth  transition  from  no  smoothing  through  1  cm  s’1  to  the 
1.5  cm  s'1  curve.  Further  smoothing  does  not  significantly  change  K(Af)  from  the  curve 
shown  for  1.5  cm  s’1.  In  what  follows  the  theoretical  results  obtained  using  a  smoothing 
equal  to  1  cm  s  1  are  given.  This  smoothing  approximately  removes  from  the  sound- 
speed  profile  variations  at  scales  less  than  a  Fresnel  zone  size  at  AFAR  (*s  20  m)  and  so 
the  selection  of  the  1  cm  s’1  scale  is  not  entirely  arbitrary.  The  sensitivity  of  the 
intensity  results  to  the  profile  smoothing  remains  an  important  problem. 

Also  shown  in  Figure  3  is  the  approximation  to  K(Al)  given  by  (3).  The  decorrela¬ 
tion  scale,  u1,  is  obtained  from  equation  (44)  of  our  companion  paper.8  The  data  appear 
to  follow  the  Gaussian  form  of  (3),  with  perhaps  a  slightly  higher  value  of  v',  at  least  as 
well  as  the  form  of  K(At)  from  more  elaborate  calculations.  Since  there  are  environ¬ 
mental  uncertainties  in  the  calculation  of  u we  feel  that  we  cannot,  from  the 
correlation-function  data,  distinguish  between  K(Al)  and  its  Gaussian  approximation. 
Another  approximation  for  K  (A/);  that  is,  |  Q  |2  from  (40)  of  our  companion  paper,  is 
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shown  as  well.  If  the  prediction  for  K(A/)  in  full  saturation  is  used 

K(A«)  =exp[— (t/At)2]  (9) 

there  would  be  a  substantial  disagreement  with  data,  since  v  is  substantially  larger 
than  i/1  (see  Table  1).  In  Figure  4,  the  observed  correlation  functions,  r(Af),  are  plot¬ 
ted  for  three  acoustic  frequencies,  with  90%  confidence  intervals  for  a  few  lags.  The 
observed  intensity  decorrelation  times  are  0.56,  0.13,  and  0.03  h  for  410,  1010,  and 
4670-Hz  respectively.  Overplotted  as  solid  curves  are  the  K(A<)  predictions  (1  cm/s 
smoothing).  At  the  90%  confidence  level  the  theoretical  curves  are  in  agreement  with 
the  observations  over  much  of  the  range  in  correlation  and  the  comparisons  are  in  gen¬ 
eral  quite  good.  The  K(A/)  curves  at  different  frequencies  are  related  by  a  simple  time 
scaling  that  is  proportional  to  <r(ln<t>)'4  with  a  the  acoustic  frequency.  The  predicted 
behavior  of  K(Af)  with  both  At  and  acoustic  frequency  is  as  observed. 

The  power  spectrum  and  the  correlation  function  are  Fourier  transform  pairs  and 
we  now  turn  to  consider  the  observations  in  the  frequency  domain.  Note  that  r(Acr)  is 
not  a  transform  pair  with  the  curves  in  Figures  4;  it  is  the  correlation  of  intensity  as  a 
function  of  frequency  lag  and  is  treated  later. 

The  power  spectra  of  the  observed  35-km  intensity  time  series  P/(w)  are  plotted  in 
Figure  5.  Standard  FFT-techniques  were  used.22  These  spectra  were  obtained  by  nor¬ 
malizing  the  series  to  a  mean  intensity  of  unity  and  then  averaging  spectra  from  50% 
overlapping  23.12-h  segments.  A  missing  5.6-h  segment  in  the  4670-Hz  data  was  filled 
using  linear  interpolation.  Each  segment  was  mean-zeroed  and  filtered  with  a  cosine 
taper  to  reduce  leakage.  The  spectra  obtained  from  the  segment  averages  were  bin 
averaged  over  0.05  logarithmic  intervals  in  frequency  w  (to  avoid  confusion  with  acous¬ 
tic  frequency  <7,  the  temporal  transform  variable  u  is  used  in  the  text);  90%  confidence 
intervals  are  shown.  The  observed  spectra  have  been  multiplied  by  2/ <  I2  >so  that  a 
comparison  can  be  made  with  the  Fourier  transform  of  K(At).  (The  integrals  of  the 
spectra  are  equal  to  two.) 

The  observed  spectra  have  an  uP  —<J~l  dependence  at  low  oj.  At  a  higher  w  (which 
is  dependent  on  acoustic  frequency)  the  spectra  begin  to  roll  off.  The  acoustic-frequency 
dependent  break  in  slope  in  the  power  spectrum  has  also  been  observed  elsewhere.11 
There  appear  to  be  possible  influences  from  the  semidiurnal  tide  (  u  =  .08  cph  )  appear¬ 
ing  as  a  bulge  at  low  u  in  the  4670-Hz  intensity  spectrum.  The  spectra  in  Figure  5  were 
calculated  in  a  way  meant  to  remove  the  longer  period  variability  and  so  low-frequency 
lines  are  smeared.  At  large  w,  the  noise  floors  are  four  to  five  decades  below  the  energy 
containing  parts  of  the  spectra. 
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The  predictions  obtained  from  the  Fourier  transform  of  K(Af)  are  overplotted  in 
Figure  5.  The  behavior  in  u  and  acoustic-frequency  a  is  very  similar  to  the  observations. 
The  measured  spectra  fall  below  the  predicted  curves  at  high  u.  Because  the  high-w 
tail  is  sensitive  to  the  smoothing  parameter  this  is  not  felt  to  be  a  failure  in  the  theory. 

To  see  the  distribution  of  spectral  energy  with  u  consider  the  variance  preserving 
plots  shown  in  Figure  6.  Equal  areas  under  the  curves  contain  equal  contributions  to  the 
integral  of  the  spectrum.  The  theoretical  curves  match  the  observations  except  very 
near  the  maxima.  They  are  within  a  factor  of  two  of  the  observations  of  the  distribu¬ 
tion  of  the  fourth  moment  in  u.  The  peak  in  the  410  Hz  spectrum  at  0.3  cph  is  not 
significantly  different  from  the  prediction  at  the  90%  confidence  level.  Finally,  we 
remark  that  the  theoretical  predictions  for  Figures  5  and  6  would  not  be  significantly 
different  if  the  Gaussian  approximation  (3)  for  K(Ai)  were  used,  if  slight  adjustments  in 
the  values  of  u  1  were  allowed. 


The  intensity  correlation  in  frequency  r(A/7),  given  by  (2),  is  plotted  in  Figure  7. 
The  short  1.25-ms,  3199-Hz  pulses  were  used  to  estimate  r(Ao-).  Each  pulse  was  Fourier 
transformed  and  it’s  spectrum  calculated  and  then  normalized  by  the  average  of  the 
pulse  spectra  observed  over  the  entire  187  h  experiment.  A  correlation  function  in  fre¬ 
quency  was  then  obtained  from  each  normalized  pulse  spectrum  and  averaged  over  all 
pulses.  The  result  is  plotted  in  Figure  7  with  the  prediction  obtained  from  |Q(Ao-)|2 
The  comparison  is  quite  good  down  to  a  level  less  than  0.2  in  r(A <r).  Also  shown  is  the 
Gaussian  approximation  to  |  Q  |2  from  (4). 

Of  interest  is  another  fourth  moment  that,  like  r(Aa),  uses  the  theoretical  quan¬ 
tity  |Q(Acr)|2  for  prediction;  it  is  the  intensity-lagged  correlation  given  by 

QO  00 

f(T)  =<fl(t  +  T)l(t)dt>/<fl*(t)dt>  (10) 

-oo  -oo 

The  lagged  correlation  coefficient  of  intensity  is  formed  for  each  1.25-ms  pulse  at  3200 
Hz.  The  average  /(r)  is  then  obtained  by  averaging  the  correlations  from  all  pulses 
and  normalizing  by  the  value  at  lag  r  =0.  Note  that  unlike  I(r),8  /(r)  does  not  depend 
upon  the  pulse  arrival  time  and  so  does  not  contain  the  wander  of  the  pulse  in  time. 

The  prediction  for  /(r)  requires  knowledge  of  the  transmitted  pulse;  equation  (15) 
of  Worcester,  Williams  and  Flatte'9  is  used. 

OO 

/(r)  f  /  /  dtedcrdo'  P(A<r.cr,</)[  |Q(A^)  |2  +  |Q(cr (II) 

-Tt  -oo 

where 
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s 


P(AwO -*(  <T + )  x-( »  -  )  *(  y  +  ^ )  -r(  ✓  -  ) 


(12) 


(~A> 


\  2  ' 

The  transmitted  pulse  shape  was  a  Gaussian  given  by  x(  t  )  =  e  20  so  that 


.Y(  a  )  oc  e  2 


(13) 


where  a  =  0.53  ms. 

Substituting  (13)  into  (12)  and  (11)  gives 


/(')-* 


r2 

-2lrL  4.p  2a2 

5(0)  +e 


(14) 


where 


g(r)  =  J  dti  |  Q  (ti)  |2  exp 


|  — %  a2  t/2+  i  w  r  j 


(15) 


The  first  term  in  (14)  depends  on  the  medium-induced  pulse  spreading,  represented 
by  the  Q  function.  The  second  term  is  a  replica  of  the  transmitted  pulse.  Note  that  at 
r=0  the  two  terms  contribute  equal  amounts.  We  may  get  some  idea  of  the  behavior 
of  /(r)  by  using  the  approximation  for  |Q|2  in  (4);  this  yields 


/(')=* 


exp 


2a2+4r,? 


exp 


2or 


(16) 


Thus,  if  Tt«a,  /(r)  is  a  replica  of  the  transmitted  pulse.  If  T0»a,  then  f(r)  has  a 
replica  of  the  transmitted  pulse  sitting  on  top  of  a  much  wider  Gaussian  with  rms 
spread  of  1.4  r0. 

However,  (3)  is  not  a  good  approximation  for  the  tail  of  /(r),  because  that  tail  is 
dominated  by  the  low-lying  eigenvalues  discussed  in  reference  8.  Thus  ^(r)  must  be  cal¬ 
culated  numerically  for  AFAR.  This  has  been  done,  and  the  result  is  shown  in  Figure  8. 
Thus  intensity  measurements  can  be  used  to  measure  the  effect  of  these  low-lying  eigen¬ 
values,  and  |Q|2  cannot  be  characterized  by  r0  only. 
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m.  PROBABILITY  DISTRIBUTIONS  AND  MOMENTS 


For  the  2.8-km  data  series  the  expectation  is  that  the  log-intensities  are  normally 
distributed.  The  time  series  of  log-intensity  were  divided  into  35-minute  segments  with 
each  segment  normalized  to  have  a  mean  intensity  of  unity.  The  log-intensities  from 
each  segment  were  then  binned  to  form  histograms.  From  an  average  of  the  histograms 
the  probability  densities  were  calculated  for  each  acoustic  frequency  with  the  results 
plotted  in  Figure  9.  This  technique  removed  dominant  oscillations  at  low  frequency  in 
the  7  hour  data  set  at  750  Hz.  There  is  very  little  difference  in  the  1500  and  2800  Hz 
distributions  whether  short  or  long  segments  are  used.  The  confidence  intervals  are  cal¬ 
culated  using  the  technique  given  in  Appendix  B  of  Spiesberger  and  Worcester.14  In  gen¬ 
eral,  the  2.8  km  log-intensities  are  normally  distributed  as  expected.  A  normal  distribu¬ 
tion  is  overplotted  for  each  series  with  the  mean  and  variance  given  by  the  observed 
distributions. 

The  observed  probability  density  functions  for  three  frequencies  from  the  35-km 
range,  are  shown  in  Figure  10.  The  axes  are  log-linear  so  that  an  exponential  would  be 
a  straight  line.  Below  the  10-3  level  in  probability  the  90%  confidence  intervals  are 
quite  large  due  to  the  small  number  of  observations.  Note  that  the  confidence  intervals 
at  small  intensity  for  410  Hz  are  broader  than  for  4670  because  of  the  longer  correla¬ 
tion  times  in  the  410  Hz  series.  The  data  fall  nearly  linearly  in  probability  space  for 
almost  two  decades  and  then  lead  into  a  tail  that  becomes  more  prominent  with 
increasing  acoustic  frequency.  This  is  the  behavior  modeled  by  the  K-distribution,"0 
shown  as  solid  curves  in  Figure  10.  The  parameter  of  the  distribution  is  the  normalized 
intensity  variance,  which  is  obtained  from  the  <I2  >  values  given  in  Table  2.  In  all 
cases,  the  data  are  well  represented  by  a  K-distribution  with  the  observed  intensity 
variance. 

The  observed  intensity  variances  are  1.24,  1.21,  and  1.47  for  410,  1010  and  4670 
Hz.  The  4670-Hz  scintillation  is  higher  than  expected.  The  values  predicted  for  7,  0.15, 
0.11  and  0.064  give  intensity  variances  of  1.3,  1.22,  and  1.13  for  410.  1010,  and  4670  Hz 
respectively.  Agreement  between  theory  and  experiment  is  reasonable  at  410  and  1010 
Hz.  The  cause  of  the  disagreement  at  4670  Hz  is  unknown,  although  it  is  known  that  7 
is  sensitive  to  the  shapes  of  the  profiles  in  sound  speed  and  buoyancy  frequency. 

The  observed  normalized  moments  of  intensity,  <I'V  >/<I  >v,  for  N  =  2  to  6  are 
given  in  Table  2  for  410,  1010,  and  4670  Hz.  The  observed  moments  increase  faster  with 
,'V  than  /V!,  (the  prediction  of  an  exponentially  distributed  intensity  series  in  full  satura¬ 
tion).  The  predictions  for  the  moments  obtained  from  (6)  and  (7)  are  tabulated  as  well 
as  those  obtained  using  the  moments  of  the  K-distribution  (8),  where  each  prediction 


has  been  normalized  to  have  the  correct  second  moment. 

It  is  very  difficult  to  estimate  errors  on  these  moment  determinations,  because  the 
error  on  <1^  > depends  on  the  true  value  of  <I2^  >  which  is  very  uncertain  when  N  is 
large.  The  approximate  numbers  of  independent  observations  of  intensity,  based  on 
experiment  duration  divided  by  the  intensity  correlation  time,  1  /u',  are  720,  2200,  and 
12.600  for  410  Hz,  1010  Hz,  and  4670  Hz,  respectively.  As  an  example  of  the  uncer¬ 
tainty  in  the  moment  error,  consider  <I5  >  If  the  intensity  probability  distribution 
were  exponential  the  error  in  this  moment  at  4670  Hz  would  be  18;  if  the  K-distribution 
applied,  the  error  would  be  2000;  if  the  K-distribution  applied,  but  the  dynamic  range  of 
the  experiment  were  limited  at  an  intensity  of  twelve  times  the  mean,  then  the  error 
would  be  about  200.  Thus,  within  the  roughly  estimated  errors,  all  the  one-parameter 
models  considered  fit  the  data  reasonably  well.  Extending  to  a  2-parameter  empirical 
model  such  as  that  given  by  Macaskill  and  Ewart,23  as  they  point  out,  does  not  seem 
necessary  to  fit  the  AFAR  data. 
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IV.  SUMMARY  AND  CONCLUSIONS 

We  have  compared  observations  of  acoustic  fluctuations  in  intensity  at  AFAR  with 
the  theory  of  acoustic  propagation  through  internal  waves.  The  predictions  are  sensi¬ 
tive  to  small  changes  in  the  sound  channel,  but  knowledge  of  the  sound  channel  and  the 
internal-wave  field  was  sufficiently  precise  to  achieve  good  agreement  for  the  intensity 
coherence  functions  of  time  and  of  acoustic  frequency,  as  well  as  their  corresponding 
spectra.  The  behavior  at  small  time  lag  is  governed  by  a  parameter  u1  that  can  be  cal¬ 
culated  from  a  double  integral  along  the  ray  from  source  to  receiver,  and  the  observa¬ 
tions  agree  with  the  calculated  v' .  The  intensity  coherence  of  acoustic  frequency  is 
grossly  characterized  by  a  parameter  t0  that  can  be  calculated  in  a  manner  similar  to 
v'. 

The  probability  density  of  intensity  has  been  observed  at  two  ranges,  one  in  the 
unsaturated  region  and  the  other  in  the  partially  saturated  region.  The  former  fits  a 
log-normal  distribution  and  the  latter  fits  a  K-distribution.  The  distribution  in  the  par¬ 
tially  saturated  region  is  characterized  by  the  microray  focussing  parameter  7.  Calcu¬ 
lation  of  7  from  the  internal-wave  spectrum  is  very  sensitive  to  the  sound-speed  profile, 
so  that  reliable  predictions  could  not  be  made  at  AFAR. 

We  are  grateful  for  a  number  of  useful  conversations  with  Dennis  Creamer  and 
Frank  Henyey  of  CSND,  and  for  the  hospitality  of  the  La  Jolla  Institute,  where  some  of 
this  work  was  done.  This  work  was  supported  by  the  Office  of  Naval  Research. 
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TABLE  I. 

Internal  wave  model  and  predicted  acoustic  parameters  at  AFAR 


Internal  Wave 

Model  Parameters 

13  1  m 

]• 

3 

n0B 

3  cph  km 

Acoustic  Parameters 

(1  cm  s'1  smoothing  of  the  background  profile) 


Short  Range 

<7q(Hz) 

103  A 

$ 

749  8 

2  3 

2  8 

1499  8 

1  2 

5  5 

2800  0 

0  63 

10 

Long  Range 

<70|HZ 

)  103  A 

<t> 

1 

*(h-M 

(h-1 

410 

47 

7  2 

15 

6  4 

3  9 

1010 

19 

18 

11 

16 

12 

4670 

4  2 

82. 

064 

73 

67 

TABLE  II. 

Observed  and  predicted  moments  <I‘V>/<I>V  at  35  km 


N  = 

o 

3 

4 

5 

6 

410  Hz 

arO.26 

±2.1 

Observed 

2  24 

8.5 

47.6 

339 

2780 

equation  (6) 

2.24 

8.5 

48 

380 

4000 

equation  (7) 

2  24 

8.2 

41 

264 

2000 

equation  (8) 

2.24 

8.3 

45 

335 

3200 

1010  Hz 

Observed 

2  21 

8.22 

43  5 

290 

2240 

equation  (6) 

2.21 

8  3 

46 

360 

3700 

equation  (7) 
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80 

40 
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equation  (8) 

2  21 
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44 
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4670  Hz 
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2  47 
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74  3 
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equation  (6) 

2  47 

11 

85 
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17000 

equation  (7) 

2  47 

10 

58 
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equation  (8) 
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FIG.  1.  The  A  — <&  diagram  for  AFAR.  The  points  in  the  unsaturated  regime,  A$*<1, 
are  from  the  2.8  km  experiment.  The  points  in  the  partially  saturated  regime,  A^M 
and  A#<1  are  for  the  different  acoustic  frequencies  over  the  35-km  path.  The  fully 
saturated  area  is  A4>>I. 


FIG.  2.  (a)  The  intensity  series  at  2.8  km.  The  interval  between  tick  marks  is  5  times 
the  mean  intensity.  The  series  are  for  750  Hz  (top  trace),  1500  Hz  (middle  trace),  and 
2800  Hz  (bottom  trace),  (b)  20  hours  of  intensity  data  for  three  frequencies  over  the 
35-km  path.  The  interval  between  tick  marks  is  10  times  the  mean  intensity.  The 
acoustic  frequencies  are  410  Hz  (top),  1010  Hz  (middle)  and  4671  Hz  (bottom). 


FIG.  3.  Observed  intensity  correlation  (dots,  with  occasional  error  bars  indicated)  in  the 
form  of  equation  (l)  for  1010  Hz  compared  with  K(A<)  calculated  using  two  smoothings: 
dy  =1.5  and  1.0  cm/s.  The  e  curve  is  also  shown.  The  |  Q  |2  curve  comes  from 

an  approximation  to  K  (At)  involving  the  two-frequency  second  moment.  9 


FIG.  4.  Intensity  correlation  functions  at  35  km;  the  dotted  curves  are  experimental 
data,  the  solid  curves  are  the  theoretical  predictions;  90^0  confidence  intervals  are 
shown  for  a  few  lags. 


FIG.  5.  Normalized  power  spectra  of  intensity  at  35  km  compared  with  the  prediction 
from  the  Fourier  transform  of  K(At).  Power  spectra  with  90^0  confidence  intervals  cen¬ 
tered  at  10-4  are  shown  on  the  410-Hz  graph.  The  almost  vertical  extra  line  segments 
at  high  frequency  are  the  tails  of  the  predictions  from  the  Fourier  transforms  of 


FIG.  6.  As  in  Figure  5,  plotted  so  that  equal  areas  under  the  curves  correspond  to 
equal  contributions  to  the  variance. 


FIG.  7.  The  correlation  of  intensity  as  a  function  of  frequency  difference  at  3200  Hz 
over  the  35  km  path.  The  dots  are  experimental  data  and  the  solid  curve  is  the 
thoeretical  prediction  The  dashed  curve  shows  the  prediction  obtained  by  setting 
U”  =0  :  that  is,  ignoring  the  sound  channel. 


FIG.  8.  The  correlation  of  intensity  as  a  function  of  time  lag  within  the  pulse.  The 
dots  are  experimental  data:  the  solid  curve  is  from  a  full  calculation  of  Q  (A<r),  that  is, 
(14);  the  dashed  curve  is  from  the  Gaussian  approximation  to  Q  (Atf),  that  is,  (16). 


FIG  9  Probability  densities  for  log-intensity  calculated  from  35-minute  segments  for 
the  three  frequencies  at  2  8  km  90°J  confidence  intervals  are  shown  (every  other  inter¬ 
val  for  750  Hz)  Normal  curves  with  mean  and  variance  from  the  data  are  plotted  as 
solid  lines. 


FIG  10  Probability  densities  for  normalized  intensity  at  35  km  90cc  confidence  inter¬ 
vals  are  shown  Exponential  distributions  are  shown  as  solid  straight  lines  K- 
distributions  with  Tf’s  from  the  observed  values  of  <I*>are  shown  as  solid  curves 


FIG  1  The  A  —  4>  diagram  for  AFAR  The  points  in  the  unsaturated  regime.  A^2<1. 
are  from  the  2  8  km  experiment  The  points  in  the  partially  saturated  regime.  A«fr">l 
and  Vt><l  are  for  the  ditferent  acoustic  frequencies  over  the  35- km  path.  The  fully 
saturated  area  is 
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FIG  2.  (a)  The  intensity  series  at  2.S  km.  The  interval  between  tick  marks  is  5  times 
the  mean  intensity.  The  series  are  for  750  Hz  (top  trace).  1500  Hz  (middle  trace),  and 
2S00  Hz  (bottom  trace  I.  |b)  20  hours  of  intensity  data  for  three  frequencies  over  the 
35-km  path  The  interval  between  tick  marks  is  10  times  the  mean  intensity  The 
acoustic  frequencies  are  310  Hz  (top),  1010  Hz  ( middle  I  and  1671  Hz  (bottom). 
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FIG.  3.  Observed  intensity  correlation  (dots,  with  occasional  error  bars  indicated)  in  the 
form  of  equation  (1)  for  1010  Hz  compared  with  K(Af)  calculated  using  two  smoothings: 
dy  =  1.5  and  1.0  cm/s.  The  e  curve  is  also  shown.  The  |  Q  |2  curve  comes  from 

an  approximation  to  K  (At)  involving  the  two-frequency  second  moment.  8 
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FIG.  4.  Intensity  correlation  functions  at  35  km;  the  dotted  curves  are  experimental 
data,  the  solid  curves  are  the  theoretical  predictions;  90*^0  confidence  intervals  are 
shown  for  a  few  lags. 
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FIG.  5  Normalized  power  spectra  of  intensity  at  35  km  compared  with  the  prediction 
from  the  Fourier  transform  of  KfA/).  Power  spectra  with  90^  confidence  intervals  cen¬ 
tered  at  10-4  are  shown  on  the  -UO-Hz  graph.  The  almost  vertical  extra  line  segments 
at  high  frequency  are  the  tails  of  the  predictions  from  the  Fourier  transforms  ol 
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FIG.  9.  Probability  densities  for  log-intensity  calculated  from  3.5-minute  segments  for 
the  three  frequencies  at  2.8  km.  QO^o  confidence  intervals  are  shown  (every  other  inter¬ 
val  for  750  Hz).  Normal  curves  with  mean  and  variance  from  the  data  are  plotted  as 
solid  lines. 


